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ABSTRACT 

In this paper we study various aspects of classical solutions to the affine Toda 
equations on a half-line with integrable boundary conditions. We begin by finding 
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the energy, and analysing the possible singularities of the field at the boundary. 
Using these constraints and extensive numerical investigations we classify the vacuum 

configurations and reflection factors for Ar Toda theories up to r = 5. 



September 1999 

1 Email Peter . Bowcock@durham .ac.uk M . G . Perkins@durham .ac.uk 



1 Introduction 



The study of (1 + 1) -dimensional integrable theories with boundaries is of interest, not 
only because many of the physical systems which they model have such boundaries, but 
also because of the more detailed algebraic structure which they introduce [1,2]. Besides 
the bulk ^-matrix, a reflection factor K is needed to describe the scattering of particles 
off the boundary. In a similar way to the case of the bulk S'-matrices [3] , various algebraic 
constraints allow postulates to be made for the exact forms of these reflection matrices 
[4-9]. One way these conjectures can then be checked is by using perturbation theory [10- 
12]. However, before this is possible we must first split the field into quantum fluctuations 
and a classical background field configuration. The background field configuration must 
satisfy one of the integrable boundary conditions. Unlike in the bulk, this can often only 
be achieved by a non-trivial solution which can be surprisingly tricky to find. Indeed such 
a lowest-energy background solution may not even exist since the boundary potential may 
destabilise a theory which is perfectly stable in the bulk. 

In this paper we shall look at one class of such models which have been recently studied 
- the real coupling affine Toda theories [13]. Some results on classical backgrounds and 
semi-classical scattering solutions have already been obtained in a number of special cases 
[14-17]. Here we present a more detailed study of such solutions, and an intriguingly 
complex picture emerges. 

There is one Toda theory for each affine Kac-Moody algebra g. The equation of motion 

is 

777^ r 

d^(j) + — ^n^V 3 "^ = 0. (1.1) 

Here we have used the notation that a 1 for % — 1 to r are the simple roots of g, the finite 
Lie algebra associated to g, and a = — where ip is the highest root of g (for untwisted 
algebras). The are the marks; integers such that 

r 

a + E n ^ = °- ( L2 ) 
i=i 

Also we have defined r to be the rank of g and introduced m the mass parameter, and 
(3 the coupling constant of the theory. From now on we shall arrange for m — (5 — 1 for 
convenience. For the theory defined on the whole line the bulk energy 

E = £dx^(^) 2 + ^W + ^.(^-l)) (1-3) 

is manifestly positive for real coupling constant. However we must modify this formula if 
we wish to consider the theory on the half-line, by restricting the above integral to the 
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physical region x < and adding a boundary potential term. Thus the energy for the 
boundary theory becomes 

E = dx (^-{d x <pf + i(9 t 0) 2 + p - l)j + H(0)U=o- (1.4) 

By varying this action, it can easily be shown that this corresponds to a theory with 
boundary condition 

W==^. d.5) 

It has been shown that to preserve the integrability of the theory we must take the 
potential to have the form [14, 18-20] 

V h {4>) = -2±A i M^+l* (1.6) 

i=0 V a i 

where the constants Ai must either all vanish (Neumann boundary conditions) or obey 
a number of restrictions which are given in [18]. For simply laced theories we have that 
Ai = ±1. 

An important point is that Ai can be positive as well as negative, and so whilst the 
theory in the bulk has positive energy, it is not clear whether the boundary potential term 
can destabilise the theory. To get a large negative boundary contribution to the energy 
we need that a 1 ■ <f> becomes large, and this in turn produces a positive contribution to the 
gradient terms in the energy. Which of these two competing factors wins is determined by 
the value of A^. In this paper we determine that one needs 

Ai < 1 (1.7) 

for the positive bulk energy to win out and for the theory to have energy bounded below 
Our approach to finding this result will be based on generalising the stability argument 
for the Sinh-Gordon model given in [15] which was based on a Bogomolny-like bound. We 
review this argument and explicitly extend it to A^ in section two of the paper. Then 
in section three, we find an implicit Bogomolny form of the energy for all the other Toda 
theories, and use this together with an analysis of the possible singularity structure of the 
solution at the boundary to prove the result (1.7). 

In the next section of the paper we embark on a detailed study of classical vacuum 
solutions to the Afp Toda theories. This continues and refines the analysis in terms of tau 
functions given in [16]. A general analysis of the equations, incorporating the singularity 
analysis of section three, is considered and a complete classification of vacuum solutions 
for different boundary conditions is given for A^F> up to r = 5. Some of the solutions in 
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[16] are excluded because of singularities in the physical region x < 0, whilst a number of 
new 'exceptional' solutions are found. 

In section five we consider the scattering solutions around the new vacuum solutions 
we have found. Again a complete classification of K matrices for up to r = 5 is given. 
The paper concludes with some conclusions. 

2 The energy on the half-line for the Sinh-Gordon theory and 

A { 2 ] Toda theory 

For the Sinh-Gordon model the parameters Ai appearing in the boundary potential (1.6) 
are known to be unconstrained by integrability. However, if we would like the Hamiltonian 
describing the theory on the half-line (1.4) to have energies bounded from below, we must 
have that Ai < 1 as shown in [15]. This result relied on being able to write the energy for 
the Sinh-Gordon theory in a Bogomolny form. Explicitly we have that the energy 

E = f° dx (-{d^f + \d t <P) 2 + (e^ + e"^ - 2)) - 2A x e*'^ - 2A e" t> /^ (2.1) 
J— oo \2 2 / 

can be rewritten as 

E = f ^^(9 t 0) 2 + ^(^0-v^e^ + v / 2e-^) 2 + 29 :r (e^ + e-^)) 
J— oo \2 2 / 

- (2A 1 e^ 2 - 2A e-^ 2 ) x=0 (2.2) 
= J\ dx Q(^0) 2 + \{dA - V2e^ + V2e-^f) 

- (2(A 1 - l)e*/^ - 2(A - l)e-*/^) x=0 . (2.3) 

The integral is non-negative, as is the boundary term if Ai < 1. Furthermore we can 
show that this condition is necessary for stability by taking the field to be of the form 

6 1 - De^ [ZA) 

where D is a constant which must be taken to be less than one for the solution to be 
non-singular in the region x < 0. This is the analytic continuation of the Sine-Gordon kink 
to the real coupling theory, and it satisfies the Bogomolny equation 

d x <j) - V2e* /V2 ~ + s/2e~^ = (2.5) 

so that the only contributions to the energy come from the boundary terms at x = in 
(2.3). By choosing the constant D to be close to one, we can take the field 0(0, t) to be 
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as large as we like at the boundary, and hence take the energy to be as negative as we 
like if (Ai — 1) > 0. Similarly we can show that that A < 1 is necessary for stability by 
considering the solution obtained by taking — > — in the above. 

We would like to generalise this result to Toda theories based on other algebras by 
finding a Bogomolny-like form 

E = dx Q(c^ - W t f + i(^) 2 + + V h {<j>i) (2.6) 

for the energy in these cases too. Comparing this to (2.1), we see that we can write the 
energy in this form provided that 

= V bulk = ±n l {e^-l) (2.7) 



^2 



i=0 



d x W = Yl w i d *<t>i (2-8) 
i=i 

and from the second equation we deduce that 

dW 

Wi = jrr (2-9) 
so that we must have the integrability condition 



dWj _ 

d<f>j d(pi ' 



(2-10) 



As usual, a solution of the Bogomolny equations d x <fii = Wi will automatically be 
a (static) solution to the Toda field equations which are the equations for the (1 + 0)- 
dimensional Toda molecule: 

- tf« + ^ E^a'e**'"* = 0. (2.11) 

P i=0 

The integrability of the Toda molecule can be used to show that any solution to these 
r second order equations must also satisfy r first order equations. The Toda molecule 
describe a finite (in fact r) degrees of freedom, and integrability can be taken in the strict 
Liouville sense; that is there exist r conserved quantities which Poisson-commute with each 
other. Note that now x plays the role of time in this system. By using the equation of 
motion (2.11) to eliminate higher derivatives, the conserved quantities which we label Hi 
can be written as functions of (pi and the momenta pi = d x <f>i. Thus we have 

i?i(&,Pi) = 7i- (2-12) 
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We shall be interested in solutions for which 4>i,pi — > as the 'time' x — * 
these are the solutions which will have finite energy in the (1 + 1) Toda system. For such 
solutions 7j = if, (0,0). For instance one of the conserved charges is the energy for the 
(1 + 0) system 

Hi = E \v\ ~ V bulk = ± \ V \ -£n, t (e^ - 1) = 0. (2.13) 

i=l Z i=l z i=0 

Together with the other r — 1 equations from conserved charges, we have r equations which 
can be used to solve for p { in terms of the We suggestively write these solutions as 

p t = d x <p t = Wi{<j>). (2.14) 

For these equations to satisfy the criteria to be Bogomolny equations for the Toda system 
we need to show that the conditions (2.8) and (2.10) hold. The first condition follows 
immediately from (2.14) and (2.13). The second condition follows from a result in classical 
mechanics that any equations derived from Poisson commuting quantities must also be 
Poisson commuting (see for instance [21]). Thus we must have that 

This procedure can be explicitly performed in the case of the affine Toda theory based 
on the algebra ■ For this case it is convenient to use variables 

m = e x ^ (2.16) 
Tii = \i-d x <f) (2.17) 

where as usual \ are the highest weights of the two fundamental representations of A 2 (i.e. 
the 3 and the 3). In terms of these variables the two conserved quantities can be written 

Hl = nl + nl -7r l7 r 2 -^-^- — + 3 = (2.18) 

U 2 Ui U\U 2 

9 9 «9 it? ir 2 fl , 

H 2 = Tl^l - Tx\lX 2 - ^7T! + ^7T 2 + = 0. 2.19 

Ui U 2 UiU 2 U1U2 

These equations can be easily solved for 7Ti and 7r 2 as follows 



/H 




- u 2 





TTi = («i-l)W (2.20) 

V uiu 2 



/1 + M1 + M2 , 001 , 

vr 2 = {u 2 -l)\ • (2.21 

V uiu 2 
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Other solutions exist but these do not correspond to real 0. From this we see that 



W = 2 P + UI+U2) \ (2.22) 

V UiU 2 

From the Bogomolny argument it now follows that the energy of any field configuration on 
the half-line is bounded by 

E > [Wt^ + V^o (2.23) 



-6^3 + 2 /(1±^±^1! -2A^- 2A 2 ^ - 2A -L= (2.24) 



w (1 + Ml + U 2 ) 3 - Ju\ - Ju\ - 1 

= -6V3 + 2-Y-: — =4 

JU1U2 

- 2(A 1 - l)-£= - 2{A 2 - 1)-^= - 2(A - 1)-^=. (2.25) 

For real 0, tij > and it is easy to show that the second term is bounded below in this 
region. On the other hand, the last terms are also non-negative for 

At < 1. (2.26) 

These conditions ensure classical stability, and indeed it can be shown that these conditions 
are also necessary. 



3 Bogomolny equations and stability for other afline Toda 

theories on the half-line 

Let us now proceed to the cases of other algebras. In principle we could follow the same 
steps as for A^ and A^ but in practice we are unable to invert the conserved quantities 
to find explicit relations between for the momenta in terms of the fields 0. However, we 
can use our knowledge of (analytically continued) static solutions to the Toda equations 
to circumvent this difficulty Any static solution must necessarily obey the conservation 
laws for the Toda molecule, and so in turn must satisfy the Bogomolny equations (2.14), 
and saturate the Bogomolny bound. It follows that the energy density for such a solution 
can be written is a total derivative of some function W. But it has been known for some 
time in the context of imaginary coupling Toda theories that this is indeed the case, and 
the explicit formula for W is given in terms of tau functions as 

W = -2±y^ (3.1) 

i=0 "i '1 
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where 



± 2 -^Hn). (3.2) 

i=0 a i 



Some clarifying remarks about the above formulae are in order. First, note that the 
equation (3.2) cannot completely specify the r + 1 tau functions, since there are only r 
components of 0; indeed we can send 

r, t ^n(f(x,t)r (3.3) 

where 

E^ = 0, (3.4) 

i=0 a i 

without affecting (3.2). This ambiguity is partially removed by fixing a particular form of 
the resulting equations of motion for the tau functions; 



7~i TV To 



2 r»T l + T?=[l[^- K ^T?)^t (3.5) 
Vi=o 



where = 2ai ■ aj/a^ is the extended Cartan matrix. We are still free to do a trans- 
formation of the form (3.3) for any / that satisfies <9 M <9^ ln(/) = 0. For time independent 
solutions this implies that 

/ = e ax+b (3.6) 

which only adds an unimportant constant to W. We complete the 'gauge fixing' by de- 
manding that t — > 1 as x — > — oo. With this choice of tau functions all the static solutions 
to the real Toda equations can be described in terms of r parameters, one real parameter 
for each node on the (non-extended) Dynkin diagram corresponding to a real fundamental 
representation, and one complex parameter for each pair of complex conjugate nodes. One 
can now consider the map between the parameters in the tau functions and the value of 
the corresponding static soliton r-component field on the boundary which is obtained by 
putting x = in (3.2). We shall simply assume that imposing the condition that the field 
is free of singularities in the physical region and and tends to zero as x — > — oo renders this 
map invertible so that we can write the parameters in terms of the values of the fields of 
the corresponding solitons at x — 0. In this way we can view W as a function of <fi as we 
did for and With this in mind, we can write the energy on the half-line of any 

field configuration in the form (2.6) so that the energy is bounded below by 

E > E bound = W(<f>)\ x=0 - W{4>)\ x =-oo + H(0)U=o (3.7) 

i.e. the energy of the static soliton configuration with a value of that coincides with that 
of our arbitrary solution at x — 0, — oo at some point in time. 
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To prove stability of Toda theory on the half-line, we must show that the right hand 
side of (3.7) is bounded from below. This can be written 



2 t[ , 2n, 



where = 2a.i ■ aj/a 2 is the extended Cartan matrix, and we have used that W((p)\ x= ^ 00 
vanishes since rl(x) — > as x — > — oo. If the energy is not bounded below, we can tune the 
parameters defining the tau functions to make E^ound arbitrarily negative. Naively this can 
occur two ways: either one of the tau functions become very large, or else becomes zero. 
In fact the first possibility never arises, since if one tunes the parameters to make the tau 
function large, the quotients appearing in E^ound tend to a finite limit. So if the energy is 
unbounded we must have a soliton solution with one or more tau function vanishing at the 
boundary and for which E^ound — —oo. 

Our approach therefore is to consider the function E^ound for such soliton solutions with 
singularities at the boundary, and determine whether the residue of the pole in x is positive 
(corresponding to infinite positive energy or negative). If the residue is always positive (or 
zero) then E^ound must be bounded below. To analyse the behaviour near the origin take 

n = ai x yi + 0(x yi+1 ) (3.9) 

and we know that Hi > and a, ^ for all i. Let S C 0, 1, ...r be the set of i for which 
Hi > 0, i.e. for which the corresponding tau function Tj vanishes at x — 0. For % e S the 
equation of motion (3.5) contains a term of the form x 2yi ~ 2 . Comparing coefficients we 
find that for 

Vifla?« = (3.10) 

j=0 A 

jZK l3 y 3 = 2. (3.11) 

j=0 

Note that the second equation implies that not all the tau functions can vanish simultane- 
ously, or in other words S ^ {0, 1, ...r} since then we could write 

o = E «i— V 1 ^' = ? niKijVj = 2 H n i = 2h ( 3 - 12 ) 

which is clearly not true. Also from both equations we see that for i e S, then for x close 
to 0, 

^ = Ut -«, „ n ( a,^)^ = (3.13) 

j j=0 /ijUjX 
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Thus the poles in Eb oun d at x = have the form 



^bound 




~ 2—i I 2 ~ i 2 

x ies \ a i a i 



2 V (2 yi 



) 



+ 0(1) 



(3.14) 



where we have used that e ai '^ > and x < to identify the correct sign for the last term. 
Thus we see that E bound will be bounded below if and only if 



(3.11) implies we can always arrange yi — 1 by arranging that only Tj vanishes. Note that 
this result implies that all simply-laced Toda theories with integrable boundary conditions 
are stable, although those with Ai = 1 are only marginally so. 

4 Vacuum soliton solutions to the Toda field equations 



Toda equations which satisfy integrable boundary equations. A large class of solutions 
was found in [16] and subsequently discussed in [17]. These solutions typically involved 
pairs of analytically continued solitons [22]. However, the analysis could not determine 
whether solutions (beyond those containing a single pair of solitons) had any singularities 
in the physical region x < 0, which would render them unacceptable since such solutions 
have infinite energy. Nor was it clear whether there were other solutions to the equations, 
potentially with lower energy and therefore the 'true' vacuum solutions. Here we carry the 
analysis a little further, utilising the singularity analysis of the previous section. Extensive 
numerical analysis indicates that generally the vacuum solutions of [16] with more than 
one pair of constituent solitons do develop singularities in the physical region and so are 
unacceptable. Moreover, we shall see that even for critical values of the soliton parameters 
where it is possible to place such singularities on the boundary, we cannot obtain acceptable 
vacuum solutions. However a number of exceptional solutions are found which provide 
acceptable vacuum configurations for A^ up to r = 5. 

Let us recall some of the notation of [16] . In the language of the tau-functions introduced 
in the last section, the boundary conditions can be conveniently written (specialising to 
the A^ series of affine Toda theories where rij = 1 and af = 2) 




In this section we elaborate on the results in [15, 16] looking for static solutions to the A^ 



+ A t e aK<p/2 = C 



(4.1) 



where 




(4.2) 
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Using this definition of C, the equations of motion imply, at the boundary x — 0, 

n - — - 7? + 2Cr; - (C 2 - l)n = (4.3) 

which simplifies in the static case to 

t>'-2Ct: + (C 2 -1)t 1 \ x=0 = 0. (4.4) 

However this equation only contains the squares of the coefficients Ai and so we no longer 
know which boundary conditions a given solution of this will obey. We will find that it is 
possible to recover this information later. 

The parameter C turns out to be proportional to the energy of the field by 

E = -2Y j n i C. (4.5) 

i=0 

It is now desirable to try to find the lowest energy (i.e. highest C) static background 
solutions to these equations. These will depend, in general, on the boundary conditions 
imposed. It is known that there are solutions to the equations which consist of solitonic 
solutions analytically continued from the imaginary-coupling theory [22, 16]. The tau func- 
tions for an iV-soliton solution can be written as: 

1 1 / jv _ \ 

Tj (x,t) =£.-£ exp £/i p ($ p + Ino/w") + £ Wq \nA^ . (4.6) 
mi=o n N =o \p=i i<p<q<N J 

Here, $ p = cr p (x — v p i) + £ p . We have the mass-shell condition cr 2 (l — f 2 ) = m 2 p , which in 
the static case implies o = m ap since we need the asymptotic form — > 1 as x — > — oo. For 
the masses m ap = 2sin(^.) and uo =exp(^-). Later it is useful to take £ p = \n(d) + ix 
where d is referred to as the position of the soliton. 

The interactions between solitons are given by the interaction constants 

A{ap a q ) = K ~ ° q ? ~ {<r P v P ~ ° q v q ) 2 ~ 4 sin 2 ^(a p - a q ) 
(<t p + °q) 2 - (vpVp + (y q v q ) 2 - 4 sin 2 ^(a p + a q )' 

It is worth noting that when we interact one soliton with another of the same type and 
velocity, we find that the interaction constant A ( - apaq ^ vanishes. Hence in the static case it 
is only possible to have at most one of each type of soliton in a background solution. 
It is convenient to make the change of basis 



k 
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which splits the tau functions up into 'charge sectors' T^. This gives us a similar equation 
to before but now for T k : 

T'l - 2CT' k + (C 2 - l)T k = (4.9) 

at x — 0. We now have a simpler set of equations: a soliton of type a resides in T a whilst 
it's interaction with another soliton b resides in T a+ 5 and so on. 

It is required that the function be real everywhere on the interval (— oo,0]. Thus 
the tau functions must be real and non-negative on this interval. Reality imposes the 
condition that each soliton must appear with its conjugate, i.e. type a with type r + 1 — a. 
The exception to this is where we have a middle soliton. The middle soliton has a = 
when r is odd, and hence if we take x = 0, vr we obtain a real tau function. However the 
restriction imposed by insisting that the tau function be non-negative is more subtle and 
will be considered case by case below. Let us look at possible background configurations 
in turn. 



4.1 = 

Looking at the boundary conditions (4.1) and remembering (1.2), it is clear that the 
solution = is only valid when 

= ]T^a\ (4.10) 

i=0 

This is clearly only true when all of the Aj, are the same sign. In addition the value of C 

is found to be ±1 for the + + ..+ and ..— boundary conditions respectively. We use 

the notation for the boundary conditions where the signs of the parameters Ai are listed 
in turn. 



4.2 Single Middle Soliton 

In the case r = odd, there is a self-conjugate soliton associated with the central spot of the 
Dynkin diagram, i.e. of type a = If we consider a background configuration consisting 
only of such a static soliton, we find; 

Tj = l + (-1)3 de 2 * (4.11) 

where we have taken x — for simplicity (x = 71 merely swaps the role of the even and 
odd tau functions). It is apparent from this equation that we must take d < 1 in order 
that none of the tau functions are negative in the physical region. 
Since 

eQ i.* = Zhljil, (4.12) 
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the boundary conditions imply that 

1 - {-lyde 



l + (-l)^e 2 - ~ Aj 



-lV'de 



2.r 



1 + (-l)J'cfe 2a; 



(4.13) 



Note that here we have used the fact that the energy of such a solution is given by C = 
m — 1 = 1. For d < 1, the modulus sign is irrelevant and we obtain the result that all 
the Ai must be positive. In the special case d — 1, all the odd tau functions vanish at the 
boundary. Consider a vanishing tau function, Tj. Since we can write the j — 1 and j + 1 
boundary conditions as 

£7-^1 = ^ — ^7^ (4.14) 

T j-i 

and 

C-^i J+1 l^ (4-15) 

respectively, then it is clear that in this case we can take Aj-i and Aj + i to be of arbitrary 
sign. This gives us a variety of consistent boundary conditions. It is worth noting also 
that since rj < and Tj > for x near to but below 0, then the boundary condition 
^ e ai-<p/2 _ q _ _j_ > q £ Qr x _^ q an j go j£ T . vgjughgg a t x = then = 1. Hence when 

d = 1, we can fit boundary conditions of the form a + a + a + a + .. where "a" denotes 
an arbitrary choice of sign. Note that taking instead x = n allows us to fit the boundary 
conditions of the form + a + a + a + ... 

4.3 TWO SOLITON SOLUTIONS 

The single middle soliton was something of a special case: let us now consider the case 
where we have a soliton, type a, and its conjugate, type a = r + 1 — a. Then the tau 
functions become; 

n = l + 2d cos ( X + e maX + A^d 2 e 2maX . (4.16) 



r + l) 

In this case it can readily be shown (using the equation for the highest occupied charge 
sector) that there are two solutions C± = m a ± 1. One can show that the lower energy 
solution C + has singularities in the physical region so we require the parity inverse solution 
C_ [16]. There is also a specific value of d corresponding to this solution (determined by 
the zero charge sector) which is given by d — 2 2 m . 

Working with the boundary conditions (4.1) as before it can be shown that the coeffi- 
cients Ai are given by 

a ■ ( (X , n(i + l)a\ (x , vr(i-l)a\\ 

= -*S* cos - + cos - + -^-^ . (4.17) 
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Note that a shift of the parameter \ by 2n/(r + 1) merely cyclically permutes the boundary 
conditions obeyed. However, by choosing x so that cos(| + ^) = we can ensure that 
Tj — > as x — > 0. So as before this allows us to take the signs of A,-_i and to be 
arbitrary. Therefore there are a number of boundary conditions consistent with a two- 
soliton solution which is singular at the boundary. 

4.4 MULTI-SOLITON SOLUTIONS 

Let us first consider the case developed in [16] where the highest charge sector occupied 
Qmax — S p a p < We shall refer to this case as having no overlapping of the charge- 
sectors, i.e. the interaction term involving all the solitons a p always resides alone in the 
Qmax charge sector. This leads to the expressions for C and d p : 

N 

c ± = ]Tm ap ±i ( 4 - 18 ) 

p=0 

d TT ™«r+ m «P 2 (4 19) 

r A _4 \m ar -m ap \ 2Tm ap 

excepting that when r is odd we are only allowed the C_ solution. 

Extensive numerical work with 3, 4, and 5 soliton solutions suggests that, under this 
regime, the only solutions which are non-zero in the region (— oo, 0) occur when a maximal 
number of tau functions vanish at the boundary x — 0. Consider a solution in which two or 
more consecutive tau functions vanish at the origin. In the language of section three, this 
means that the S contains at least two neighbouring nodes. Let us rewrite the conditions 
(3.10) and (3.11) for the special case of theory. These give 

2y i -y i+1 -y i _ 1 = 2 (4.20) 
^± = V, (4.21) 

a i 

Now in the case where S contains at least two consecutive nodes it follows immediately 
from (4.20) that y^, the order of zero of r« at x = must be greater than one. However, 
using (4.20) we see that for i e S, both sides of the boundary condition (4.1) 

(4.22) 

go as 1/x as x — > 0. Hence, since C (the energy) is finite, the coefficients of these leading 
order terms must match. This requires, using (4.21), 

Vi = AVVi ( 4 -23) 
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which cannot be satisfied for y, L > 1. This tells us that only non-consecutive tau functions 
tending to zero at the boundary are allowed. This result is important since numerical in- 
vestigation suggests that, under the non-overlapping charge sector regime, it is not possible 
to find multi-soliton solutions which are both regular in the physical region and obey this 
restriction. However, a proof of this conjecture is lacking. 

4.4.1 Overlapping charge sectors, Q max > ^ 

It is however sometimes possible to obtain acceptable multi-soliton solutions in cases where 
the charge sectors overlap. The first case where this is possible is r = 4. Consider a solution 
containing all the possible solitons of the theory. Then look at the charge sectors: 

Charge sector Soliton combinations 
-2 2 211 12 

-1 T 122 12 

IT 22 1T22 

1 1 122 T2 

2 2 21TT2 

Now we must apply the usual charge sector equations (4.9). In the non-overlapping 
case, the highest occupied charge sector yields an equation solely in C. Here, however, we 
obtain an equation in d\ and C which must be solved simultaneously with the other charge 
sector equations. The three equations from (4.9) are: 

C 2 -l + d\A lT [{C - 2m0 2 - 1] + d\A 22 \{C - 2m 2 ) 2 - 1] 

+ ^^(A 12 A 12 ) 2 A lT A 22 [(C - 2(mi + m 2 )) 2 - 1] = (4.24) 

(C- mi ) 2 -l + d 2 2 A 12 A 12 A 22 [(C - (mi + 2m 2 )) 2 - 1] 

+ d 2 e l{x2 ~ 2xi) A 12 [(C - (mi + m 2 )) 2 - 1] = (4.25) 

and 

(C-m 2 ) 2 -l + <P 1 A 12 A 1 *A lJ [(C-(m 2 + 2m 1 )) 2 -l] 

+ d ie - i{xi+2x2) A 12 [{C - (mi + m 2 )) 2 - 1] = 0. (4.26) 

If we take e^ X2 ~ 2xi ^ and e~^ Xl+2x2 ^ to be real we can find another real solution to 
these equations. We shall call this the "exceptional" solution. This has energy C = 
mi + m 2 — \/5 and we can once again determine the boundary conditions which it obeys. 
What is remarkable is that in this case, this exceptional solution is regular in the physical 
region, has non-consecutive zeroes of the tau functions at x — 0, and fits exactly those 
boundary conditions not covered by the non-exceptional solutions. 

We can also use this technique for the case r = 5 where, despite greater difficulties with 
finding solutions to the four simultaneous equations, we obtain similar results. 
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4.5 Vacuum Solution Results 

Before presenting our results for the static background solutions for the series of affine 
Toda field theories up to r = 5 below, let us quickly recap the allowed background solutions. 
Of the non-exceptional solutions, there are few allowed background configurations. In the 
case r = even, the only allowed background solutions contain either no solitons or a 
soliton/conjugate soliton pair. When r = odd we are also allowed the case of a single 
middle soliton. In any of these cases, the singularity may reside at the boundary x = 
0, allowing additional choices in the boundary conditions obeyed. However, when r > 
4, these configurations do not span all the possible boundary conditions. To find the 
vacuum solutions for these remaining boundary conditions, it is necessary to consider the 
"exceptional" solutions. In the cases considered these exceptional solutions completely 
cover the boundary conditions not covered by the non-exceptional configurations. It is 
expected that similar results will extend beyond r = 5 although a proof is lacking. 

It should be noted that, in many cases, more than one possible background configuration 
fits with a particular boundary condition. When this occurs, the vacuum solution is of 
course that background configuration which has lowest energy. 

Our results up to r = 5 are tabulated in table 1. These results are ordered within each 
theory with decreasing C, or equivalently, increasing energy. 

5 Classical Scattering Solutions 

Besides classical backgrounds, more information can be obtained from the classical equa- 
tions of motion. The classical limit of the reflection matrix K a can be found by considering 
solutions which correspond to infinitesimal classical fluctuations around the background 
configuration [15]. 

Suppose the field is perturbed about the vacuum 

(f)(x,t) = (f) vac (x) +e(x,t). (5.1) 

Then by considering the equations of motion we find that, in the limit x — > — oo, e(x, t) 
may have the asymptotic form 

e(x, t) - p a e- iEt (e ipx + K a e' ipx ), (5.2) 

consisting of two superposed incoming and outgoing states. Here p a is a mass-eigenstate; 
an eigenvector of the mass matrix M 2 = ^a 1 <g) a 1 . Solutions of this form which obey 
the boundary conditions can be generated by considering non-static two soliton solutions, 
where now 

t) = a(x - vt) (5.3) 
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Boundary Condition 


Solitons 


C 


Singular at x = 


+ + + 


none 


1 


no 


- + + 

- - + 


IT 




no 

yes 




none 


-1 


no 



r = 2 



Boundary Condition 


Solitons 


C 


Singular at x = 


+ + ++ 


none 


1 


no 


- + ++ 

- + - + 


2 


1 


yes 
yes 


- - ++ 
+ 


IT 




no 

yes 




none 


-1 


no 



r = 3 



Boundary Condition 


Solitons 


C 


Singular at x = 


+ + + + + 


none 


1 


no 


- + + + + 

- + - + + 


22 




no 

yes 


- - + + + 

- - + - + 


1122 




yes 
yes 


+ + 


IT 




no 

yes 




none 


-1 


no 



r = 4 



Boundary Condition 


Solitons 


C 


Singular at x = 


++++++ 


none 


1 


no 


-+++++ 






yes 


-+-+++ 


3 


1 


yes 


- + - + -+ 






yes 


- + + -++ 






no 


-- + + ++ 






yes 


-- + + -+ 


22 


Vs-i 


yes 


-- + -++ 






yes 


-- + --+ 






yes 


+ ++ 

+ -+ 


1T223 


VS + 2- v /2(2 1 /3 + 2 2 /3 + 2) 


yes 
yes 


++ 

+ 


IT 





no 

yes 




none 


-1 


no 



r = 5 



Table 1: Constituent solitons and energy parameters C for the vacuum solutions of 
afflne Toda field theory on a half-line with integrable boundary conditions, given up to 
r = 5. 
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is defined by 

a = ±ip (5.4) 

av = iE. (5.5) 

Hence we must add these two solitons both of type a to the background configuration and 
determine K a by insisting that the resulting tau functions satisfy 

f . _ T » + 2CV; - (C 2 - l) n \ x=0 = (5.6) 

or of course equivalently from the charge sector point of view 

f t - TV + 2CT[ - {C 2 - 1)T^U =0 = 0. (5.7) 

It is then very simple to determine the values of K corresponding to different boundary 
conditions. 

Although they have already been set out in [16], let us briefly restate the results for 
the classical reflection factors in the cases which have non-exceptional vacuum solutions. 
In fact, we shall see that this is not necessarily relevant to cases which are singular at the 
boundary x = 0. Solving the equation for the highest occupied charge sector reveals that 

_ 2ip + ml " A*>jj>) 
a 2i P -mlj-} 1 A° b >(-p) { } 

where a is the particle being scattered and bj are the solitons present in the background. N 
refers to the number of soliton species (not including conjugates) present in the background. 
The only exception to this is the flat background case + + ...+, where 

K. = (5.9) 
2ip + m l a 

It has been checked in [16] that these reflection factors obey the classical reflection boot- 
strap equation. 

There is however a little more subtlety involved when we consider the cases where one 
or more of the background tau functions vanish at x — 0. In this case it is not so clear 
that a linear perturbation in — > <fi + e and in the tau functions T{ — > r« + are equivalent. 
However, with certain restrictions on the 6j (namely that if Tj — > as x — > then q — > 
as quickly or faster in order to result in a finite perturbation e) we can achieve the same 
results. Since the boundary conditions tell us that from the tau function point of view 

C -7 = *JW <"°> 
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then if we take i so that Tj+i goes to zero whilst Tj_i and do not, then clearly we 
require the term linear in e in r i+1 to go to zero at x — 0. In fact, this is almost the 
same restriction as the requirement that the perturbation e is finite. This condition gives 
us another equation which must be satisfied by the reflection factor K. It is simplest to 
consider the case where we take the parameters e %x to be real, so that the vacuum solution 
tau functions are symmetric about r (i.e. r r _ i+1 = Tj) — the other boundary conditions 
are of course related by the cyclical symmetry. Let us illustrate our point by considering the 

case of A 2 ^ for the boundary conditions H or + H — which contains two static solitons 

in the background solution. Here, let us take x — 71 so that, as far as the background 
solution is concerned, r — > as x — > and T\ = r 2 . Now adding in the perturbation, we 
require the right-hand sides of the two boundary conditions 

C-* = *Jf (5.H) 



and 



r 2 




C-^ = A 2 J^ (5.12) 



(evaluated at x — 0) to be equal up to 0(e), with a possible sign difference between A 1 
and A 2 . Hence we expect 

Al T l = A2 Il (5.13) 

up to 0(e). In actual fact, for A 1 = —A 2 , this equation is satisfied identically for all K 
and hence the same scattering data as we found before for the non-singular background 
also obeys this boundary condition. However, for A 1 = A 2 we require both sides of the 
equation to vanish. Imposing this condition leads to a somewhat unaesthetic non-diagonal 
scattering solution also found by Corrigan [23]. Diagonal scattering is expected in A^ 
as the conservation of spin even charges (which distinguish between particles and their 
conjugates) implies that a particle should be reflected back into itself [15], and is moreover 
implied by the bootstrap equations. To see this, consider the bootstrap equation depicted 
in Figure 1 for the non-diagonal case. Consider the case where c — 1. Then non-diagonal 
scattering implies that c can be reflected as a particle of type 1 or 2. This particle will 
then decompose into two particles a and b, either of type 1 or 2, but in any case a = b. 
However, the RHS implies that c will split first into two particles of type 2, which are then 
each reflected as types 1 or 2. In this case there is no requirement that a = b. Hence the 
scattering must be either diagonal (as generally believed for real coupling Toda theory) or 
completely off-diagonal (a possibility suggested by the recent paper by Delius [24]). This 
is rules out the unaesthetic non-diagonal scattering solution for the case A 1 — A 2 . 
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Hence we find that there is no consistent classical scattering data for the boundary 
condition H . 

It is not difficult to extend the above argument to include all vacuum solutions in 
affine Toda theory consisting of two static solitons. When the singularity is placed at the 
boundary it can be shown that, for any i G S, we obtain an analogous equation to (5.13), 
involving the boundary parameters A^i and A i+1 . This relation is only solved for diagonal 
scattering if we takes these two parameters to be opposite. This rules out a large subset of 
the integrable boundary conditions. Whilst it would appear that there is no problem with 
singular background solutions containing only one static middle soliton (here, the signs of 
the (say) odd i boundary parameters remain arbitrary since the 0(e) term of ^ vanishes), 
we also find that the 'exceptional' solutions of A^ and A^ have problems. For example, 
using the exceptional solution of A±~\ we find that there are two vanishing tau functions, 
Ti and r 4 (preserving the symmetry around r ). Then, as noted earlier, we require the 
linear perturbations of both of these vanishing tau functions to vanish at x — 0. But 
each of these requires a different K. A similar conclusion is reached when considering the 
A^ exceptional solution. So in these cases, it seems that no diagonal scattering solutions 
are possible which go as p a as x — > — oo. It would therefore appear that a large subset 
of the integrable boundary conditions do not admit suitable classical scattering solutions. 
The ones that do fall into one of three categories: they are consistent with a background 
configuration which either contains no solitons, a single middle soliton (in the case r = odd) 



19 



or a soliton/conjugate soliton pair where the signs of the boundary parameters on each 
side of a vanishing tau function are opposite. 



5.1 Classical Scattering Results 



We use the usual notation for expression of the reflection factors. We define 




sinh(| + 



2(r+l) 



(5.14) 



sinh(| — 



2(r+l) 



The results for the reflection factors for r up to 5 are given in table 2. 



6 Conclusions 



In this paper we have presented a number of results about classical solutions to Afline Toda 
field theory on a half-line. Extending what had previously been achieved for the special 
case of sinh-Gordon [15], we used a Bogomolny argument to ascertain which boundary 
conditions lead to stable theories with bounded energy on the half-line. This class includes 
all integrable boundary conditions for simply-laced afline Toda theories, although many lie 
on the boundary of stability. 

However, it turns out that there are still difficulties with providing a classical picture 
of scattering in some cases. This was found by considering the vacuum and scattering 
solutions for low rank affine Toda field theories in great detail in order to attempt to 
shed some light on what happens in general. We have found many unusual characteristics 
of these theories. It was found that although acceptable vacuum solutions can be found 
for all the "integrable" boundary conditions, not all of these appear to admit classical 
scattering data consistent with the reflection bootstrap equation. Another problem is that 
many of the classical background solutions contain parameters or moduli describing equally 
valid degenerate ground states. Exact solutions which correspond to excitations of these 
'zero-modes' can be found, but they are problematic in that they become singular in some 
finite time. If analytically continued beyond this time the solutions then satisfy changed 
boundary conditions. 

Finally, it would be interesting to compare these results with those of [25,26] where 
quantum reflection factors for ATFT were considered. 
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Boundary Condition 


K U K T 


+ + + 


-(1)(2) 


- + + 


(5)(f) 2 (§) 
(D(2) 


- - + 


not diagonal 




1 

(D(2) 


r = 2 



Boundary Condition 


KuK- 


^2 


+ + ++ 

- + ++ 

- + -+ 


-(1)(3) 


"(2)(2) 


- - ++ 


(2Y 
(D(3) 


(l) il (3) a 


_ _ + 


not diagonal 


not diagonal 




l 

(1X3) 


i 

(2)2 


r = 3 


Boundary Condition 


KuK- 




+ + + + + 


-(1)(4) 


-(2)(3) 


- + + + + 




(§)(§)"(*) 


(D(4) 


(2)(3) 


■ - + - + + 


not diagonal 


not diagonal 


- - + + + 

- - + - + 


not diagonal 


not diagonal 


+ + 


(I) 2 


(iX§X0(f) 


(*)(1)(4)($) 


(2)(3) 




not diagonal 


not diagonal 




l 


l 




(D(4) 


(2) (3) 


r = 4 



Boundary Condition 




K 2 ,K^ 


K 3 


++++++ 
-+++++ 
- + - + ++ 
-+-+-+ 


"(1)(5) 


-(2)(4) 


"(3) 2 


- + + -++ 
-- + + ++ 


(2)(4) 
(1)(5) 


(D(3) 2 (5) 
(2) (4) 


(2) 2 (4) 2 
(3) 2 


--++-+ 
--+-++ 
-- + --+ 


not diagonal 


not diagonal 


not diagonal 


+ ++ 

+ -+ 


not diagonal 


not diagonal 


not diagonal 


++ 


(3) a 


1 


(i)W 


(DW 


(3) 2 


-+ 


not diagonal 


not diagonal 


not diagonal 




l 


l 


l 

(3)2 




(1)(5) 


(2)(4) 



Table 2: Classical reflection factors for affine Toda field theory with integrable bound- 
ary conditions, given up to r = 5. 
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